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1 Introduction 



Infrared (soft and collinear) singularities appear in the calculation of multiparton QCD 
matrix elements. Although the singularities cancel in the evaluation of inclusive cross sec- 
tions, their factorization properties are at the basis of many important tools in perturbative 
QCD applications to hard-scattering processes [0]. 

At the leading order in the QCD coupling, as, the structure of the infrared singularities 
is well known to be universal. It is embodied in process-independent factorization formulae 
of tree-level 0-§| and one-loop |7HT0| amplitudes. These factorization formulae have played 
an essential role in the setting up of completely general algorithms |5|j7| JTT|JT^] to handle and 
cancel infrared singularities, when combining tree-level and one-loop contributions in the 
evaluation of jet cross sections at the next-to-leading order (NLO) in perturbation theory. 

The extension of these general algorithms at the next-to-next-to-leading order (NNLO) 
is at present one of the main goals to improve and precisely quantify the theoretical ac- 
curacy of perturbative QCD predictions. To this purpose we need to compute two-loop 
matrix elements [|i3] -|T5ll and to understand the structure of the infrared singularities of 
two-loop, one-loop and tree-level amplitudes at The singular behaviour of two-loop 

QCD amplitudes has been discussed in Ref. JT6|. The soft and collinear limits of one- 
loop amplitudes have been derived in Refs. [IT], P|. The soft, collinear and soft-collinear 



singularities of tree- level amplitudes have been studied in Refs. |19| , |20|j , pT|[22f and [^T 
respectively. 

The purpose of this paper is twofold. We consider tree-level matrix elements and 
present general techniques to compute their infrared singularities and to derive infrared- 
factorization formulae to any order in as- We apply these techniques to the explicit calcu- 
lation of all the relevant infrared factors at 0(a\). 

Our general method exploits the universality properties of soft and collinear emission 
and consists in directly computing process-independent Feynman subgraphs in a physical 
gauge. We use power-counting arguments p3[^4| and the eikonal approximation || to 
treat the collinear and soft limits, respectively. We show how the coherence properties of 



QCD radiation [25] can be used to deal with the mixed soft-collinear limit in terms of the 



collinear and soft factorization formulae. 

Most of the explicit results at presented in this paper were first obtained by 

Campbell and Glover [31]. The strategy followed in Ref. [^T]] was to take universal fac- 



torization for granted and thus to extract the (9(ag)-singular factors by performing the 
corresponding limits of a set of known matrix elements. We confirm their calculations by 
using a completely independent method. We also extend their results by considering the 
emission of a soft fermion pair and by fully taking into account spin (azimuthal) correlations 
in the collinear limit. The extension to azimuthal correlations is essential to apply some 
general methods to perform exact NNLO calculations of jet cross sections. For instance, the 
subtraction method [||, [12|] works by regularizing the infrared singularities of the tree-level 
matrix element by identifying and subtracting a proper local counterterm. Thus, the study 



of the azimuthally averaged collinear limit [3JJ is not sufficient for this purpose. 
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The knowledge of the infrared structure of multiparton amplitudes is also important for 
other perturbative QCD applications. The leading-logarithmic (LL) parton showers, which 
are implemented in Monte Carlo event generators U to describe the exclusive structure 
of hadronic final states, are based on the (9(as)-factorization formulae supplemented with 
'jet calculus' techniques and colour- coherence properties |3], ^5|. Analytical techniques 



to perform all-order resummation of logarithmically enhanced contributions at next-to- 
leading logarithmic (NLL) accuracy [27| rely on the factorization properties of soft and 



collinear emission. The results on infrared factorization presented in this paper can be use- 
ful to improve parton-shower algorithms and resummed calculations beyond their present 
logarithmic accuracy. 

The outline of the paper is as follows. We start in Sect. |2] by studying the collinear 
behaviour. After reviewing the known factorization formulae at O(as), we discuss the 
kinematics of the triple collinear limit. Then, in Sect. |2.3| , we present our derivation of 
factorization for the multiple collinear limit at any perturbative order. Finally, in Sect. [2.4|, 
we perform the explicit calculation of the spin-dependent splitting functions at 0(a\). Our 
results for the splitting functions were anticipated in Ref. f22|. In Sect. |] we study the soft 
behaviour. We first review the known results at 0(as) and then, in Sect. |3.2| , we compute 
the emission of a soft qq pair at Section |3]3| is devoted to double gluon emission: 

we present the corresponding soft current and obtain a compact expression for its square. 
Factorization for the mixed soft-collinear limit at and at higher perturbative orders 

is discussed in detail in Sects. T~3 and \i.fy In Sect. |] we summarize our results. In general, 



soft factorization formulae involve colour correlations. At 0(a\) these correlations cancel 
in four- and five-parton matrix elements. The explicit expressions for these particular cases 
are given in the Appendix. 



2 The collinear behaviour 



2.1 Notation and collinear factorization at O(as) 

We consider a generic scattering process involving final-state^ QCD partons (massless 
quarks and gluons) with momenta VitVii — Non-QCD partons (7*, Z°, W ± , . . .), car- 
rying a total momentum Q, are always understood. The corresponding tree-level matrix 
element is denoted by 

Ml\>%-;: sl ' S2 '--iPiiP2i---) , (1) 

where {c 1; c 2 , . . .}, {si, s 2 , . . .} and {ai, a 2 , . . .} are respectively colour, spin and flavour 
indices. The matrix element squared, summed over final-state colours and spins, will be 
denoted by |.M aija2j ...(j>i,P2, . . .)| 2 . If the sum over the spin polarizations of the parton ai 
is not carried out, we define the following 'spin-polarization tensor' 



V^ipi, •••)- E E m:,::,;;::--^- (!>..,>,. . . .) {M^^-tp^ . . .) ]t 

spins ^sijSj colours 



(2) 



§The case of incoming partons can be recovered by simply crossing the parton indices (flavours, spins 
and colours) and momenta. 
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We work in d = 4— 2e space-time dimensions and consider two helicity states for massless 
quarks and d — 2 helicity states for gluons. This defines the conventional dimensional- 
regularization (CDR) scheme of both ultraviolet and infrared [2jj divergences. Thus, 



the fermion spin indices are s = ±1, while to label the gluon spin it is convenient to use 
the corresponding Lorentz index /i = 1, . . . , d. The <i-dimensional average of the matrix 
element over the polarizations of a parton a is obtained by means of the factors 

2<^' ( 3 ) 
for a fermion, and (the gauge terms are proportional either to p^ or to p u ) 

J^2 d ^(p) = 2(1- e ) (~^ + gauge terms ) ( 4 ) 

with 

- <T <^(p) = d-2, jf d^(p) = d^(p) p» = 0, (5) 
for a gluon with on-shell momentum p. 

The singular collinear limit at 0(«s) is approached when the momenta of two partons, 
say pi and P2, become parallel. This limit can be precisely defined as follows: 



z 2p ■ n 1 — z 2p ■ n 



k 



813 = 2j?l • P2 = 7Z V , fcj. . (6) 

z(l - z) 

In Eq. (Q) the light-like (p 2 = 0) vector p^ denotes the collinear direction, while is 
an auxiliary light-like vector, which is necessary to specify the transverse component k± 
[k\ < 0) (A;_i_ • p = k± ■ n = 0) or, equivalently, how the collinear direction is approached. 
In the small-fcx limit (i.e. neglecting terms that are less singular than l/k\), the square of 
the matrix element in Eq. ([I]) fulfils the following factorization formula [|I[]: 

\M ai>a2 ,...(Pi,P2, • • Ol 2 ^ — 4vr/i 2e a s T Q ->, • • •) ^ 2 (^ k ^ e ) ' ( 7 ) 

Sl2 

where /x is the dimensional- regularization scale. The spin-polarization tensor T* s ' (p, . . .) is 
obtained by replacing the partons a\ and 02 on the right-hand side of Eq. (|2|) with a single 
parton denoted by a. This parton carries the quantum numbers of the pair ai + 02 in the 
collinear limit. In other words, its momentum is p M and its other quantum numbers (flavour, 
colour) are obtained according to the following rule: anything + gluon gives anything, and 
quark + antiquark gives gluon. 

The kernel P aia2 i n Eq. (0) is the <i-dimensional Altarelli-Parisi splitting function [0. It 
depends not only on the momentum fraction z involved in the collinear splitting a — > ai+a 2 , 
but also on the transverse momentum k± and on the helicity of the parton a in the matrix 
element Ai%'"' s '"'(p, ■ ■ ■)■ More precisely, P aia2 is i n general a matrix acting on the spin 
indices s, s' of the parton a in the spin-polarization tensor T^ s '(p, . . .). Because of these 
spin correlations, the spin-average square of the matrix element .M£;"' s '"'(p, • • •) cannot be 
simply factorized on the right-hand side of Eq. (0) . 
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The explicit expressions of P 0ia2 , for the splitting processes 

a(p) -> ai(zp + fc ± + 0(A;1)) + a 2 ((l - z)p - fc± + 0(^1)) , 
depend on the flavour of the partons ai, a 2 and are given byQ 

'1 + z 2 



P°°(z,k ± ;e) = P™(z } k ± ;e)=6 ss ,C F 



el 



(8) 



(9) 



5 S s' Cp 



1 + (1 



— 



(10) 



P 9 7(z,fc ± ;e)=P^,£; ± ;e) = T R 



P£(z,k ± ',e) = 2C A 



-9 



[IV 



+ 



1 - Z 



k>ik u 



2(l-e)z(l-z) 



where the SU(N C ) QCD colour factors are 



2N r 



C A = N C , T, 



R 



(11) 

(12) 
(13) 



and the spin indices of the parent parton a have been denoted by s, s' if a is a fermion and 
fi, v if a is a gluon. 

Note that when the parent parton is a fermion (see Eqs. @ and ([10]) ) the splitting 
function is proportional to the unity matrix in the spin indices. Thus, in the factorization 
formula (|7]), spin correlations are effective only in the case of the collinear splitting of a 
gluon. Owing to the fc^-dependence of the gluon splitting functions in Eqs. fllTf ) and flT2|), 
these spin correlations produce a non-trivial azimuthal dependence with respect to the 
directions of the other momenta in the factorized matrix element. 

Equations (|9|)-(|l2|) lead to the more familiar form of the rf-dimensional splitting func- 
tions only after average over the polarizations of the parton a. The d- dimensional average 
is obtained by means of the factors in Eqs. (|3|) and ([|). Denoting by (P aia2 ) the average of 
Paia 2 over th e polarizations of the parent parton a, we have: 



(P qg (z;e)) = e)) = C F 



1 + z 2 



e(l-z) 



(14) 



>P gq (z;e)) = (P gq (z;e)) =C F 



l + (l-z) 



— ez 



(15) 



{P qq (z;e)) = (P qq (z;e)) 



-R 



2z(l-z) 
1 - e 



(16) 



^The e dependence on the right-hand side of Eqs. (||)-((T^) refers to the CDR scheme used throughout 
the paper. A detailed discussion of the regularization-scheme dependence of the collinear splitting functions 
at O(as), including the corresponding explicit expressions, can be found in Ref. ifiof . 
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>P gg (z;e)) = 2C; 



1-z 



1 ~ z 
+ + z(l 



In the rest of this section we are mainly interested in the collinear behaviour of the tree- 
level matrix element Ai{p±, . . .) in Eq. ([!]) at At this order there are two different 



collinear limits to be considered [21]. 



The first limit is approached when two pairs of parton momenta, say {pi,pz} and 
{p3,P4}, become parallel independently. In this case collinear factorization follows from 
the straightforward iteration of Eq. (0): the ensuing factorization formula simply contains 

the product of the two splitting functions Pallt 12 and Pa^al 34 - 



In the second limit, three parton momenta can simultaneously become parallel. This 
triple collinear limit is discussed in the following subsections. 



2.2 Kinematics in the triple collinear limit 

We denote by Pi,P2 and p% the momenta of the three collinear partons. The most general 
parametrization of these collinear momenta is 

tf = Xi p» + k» ±l -^-^, i = 1,2,3, (18) 
Xi 2p ■ n 

where, as in Eq. ([]), the light-like vector p M denotes the collinear direction and the auxiliary 
light-like vector n M specifies how the collinear direction is approached (k±i-p = k±i -n — 0). 
Note that no other constraint (in particular J2i x i 7^ 1 an d J2i k±i ^ 0) is imposed on 
the longitudinal and transverse variables X{ and k±i. Thus, we can easily consider any 
(asymmetric) collinear limit at once. 

Note, however, that the triple collinear limit is invariant under longitudinal boosts along 
the direction of the total momentum pj* 23 = p± + p% + Ps- Thus, the relevant kinematical 
variables are the following boost-invariant quantities 

- " w~~r > (19) 



K = ^L-^rhr £*i • (20) 

Note that these variables automatically satisfy the constraints J2i=i z % — 1 an d k% — 0, 
so that only four of them are actually independent. 

In terms of the longitudinal and transverse variables introduced so far, the two-particle 
sub-energies are written as 

Sij = (Pi+Pj) 2 = -XiXj (— - — J = -ZiZj (— - — J . (21) 
It is also convenient to define the following variable ty^ 

/,,/, = 2 ZtS]k - ZjSlk + Slj . (22) 

Z ^ I Z j 2-"i I 2s j 
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2.3 Power counting and tree-level factorization at any order 



In the triple-collinear limit, the matrix element squared |-A4 aij a 2 ,a3,...(Pi> Pi-, P3, ■ ■ -)\ 2 has the 
singular behaviour \M ai ,a 2 ,a 3! ...(Piy P2, P3, ■ ■ -)\ 2 ~ V ( ss ')> where s and s' can be either two- 
particle (sij = (pi + Pj) 2 ) or three-particle (S123 = (pi + P2 + P3) 2 ) sub-energies. To define 
the collinear limit more precisely, we can rescale the transverse momenta k±i by an overall 
factor A: 

k ±i -> A k ±i , (23) 
and then perform the limit A — > 0. In this limit the matrix element squared behaves as 

\M aim J 2 ^0(l/\*) + ... , (24) 

where the dots stand for less singular contributions when A — > 0. We are interested in 
explicitly evaluating the dominant singular term 0(1/ X 4 ). 

To study this singular behaviour we use power-counting arguments and the universal 



factorization properties of collinear singularities. The method |24| is completely general: 
it is process-independent and can be applied to any multiple collinear limit a — > a\ . . . a m 
at the tree level (i.e. at 0(ct§~ 1 )). Thus, we shall discuss the most general case with m 
collinear partons. 

We shall show that in the multiple collinear limit a —>■ a\ . . .a m , the matrix ele- 
ment squared \M. ai ,...,a m ,...{pii ■ ■ ■ iPm, ■ ■ -)\ 2 still fulfils a factorization formula analogous 
to Eq. (0), namely 

T£(xp, . . .) P:t am , (25) 

Sl...m j 

where s\,„ m = (pi + . . . + p m ) 2 is the m-particle sub-energy and x = YhLx x %- 

As in Eq. (|^), the spin-polarization tensor TJ* S ' (xp, . . .) is obtained by replacing the 
partons ai,...,a m with a single parent parton, whose flavour a is determined by flavour 
conservation in the splitting process. More precisely, a is a quark (antiquark) if the set 
{ai, . . . , a m } contains an odd number of quarks (antiquarks), and a is a gluon otherwise. 



The factorization formula (25) takes into account all the dominant singular contributions 
in the multiple collinear limit, that is, all the contributions that have the scaling behaviour 
(l/A 2 )" 1-1 under the scale transformation in Eq. (p3|). Relative corrections of 0(A) are 
systematically neglected on the right-hand side of Eq. (|25|). 

The m-parton splitting functions P ai ...a 3 are dimensionless functions of the parton mo- 
menta pi, . . . , p m and generalize the Altarelli-Parisi splitting functions in Eq. (0). Owing to 
their invariance under longitudinal boosts along the collinear direction, the splitting func- 
tions can depend in a non-trivial way only on the sub-energy ratios Sy/si... m and on the 
longitudinal- and transverse-momentum variables Zi and fcj defined by the generalization 
of Eqs. ( |T5D and fl2H| ) to the m-parton case. 

The spin correlations produced by the collinear splitting are taken into account by the 
splitting functions in a universal way, i.e. independently of the specific matrix element on 
the right-hand side of Eq. (|25|). 
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In the case of the splitting processes that involve a fermion as parent parton, we find that 
spin correlations are absent. We can thus write the corresponding spin-dependent splitting 
function in terms of its average (P ai ...a m ) over the polarizations of the parent fermion a: 

P S ata m = 5 SS ' (Pa,..aJ • (26) 

This feature is completely analogous to the 0(a$) case and follows from helicity conserva- 
tion in the quark-gluon vector coupling. 

In the case of collinear decays of a parent gluon, however, spin correlations are highly 
non-trivial. 

Note also that the splitting functions for the collinear decay of an antiquark can be 
simply obtained by charge-conjugation invariance from those of the corresponding quark 
decay process, i.e. (P ax ...a m ) = (Pa^aJ) ■ 



The method used to derive these results exploits the basic observation |23[ that in- 
terfering Feynman diagrams obtained by squaring the amplitude M.{pi, ■ ■ ■ ,p m , ■ ■ ■) are 
collinearly suppressed when computed in a physical gauge. Thus, in the evaluation of the 
multiple collinear limit we can write 

M., (/„ p )- * [M^X Pl + ...+ Pm ,...)f V^l a J Pl ,..., Pm ) 

■ M^lipt + .-.+Pm,. ..) + ••• • (27) 

The first term on the right-hand side of Eq. (^) corresponds to the non-interfering Feyn- 
man diagrams in Fig. while the dots stand for subdominant contributions coming from 
interferences (see e.g. the diagram in Fig. ||). The superscripts (n) denote that the various 
terms are evaluated in a physical gauge. To simplify the calculation it is convenient to 
choose the axial gauge n ■ A — 0, where the gauge vector n M coincides with the auxiliary 
light-like vector used in Eq. (|18D to parametrize the collinear kinematics. The corresponding 
gluon polarization tensor dff\ is 



(») 



where q is the gluon momentum. 

The summation over spin and colour indices is understood on the right-hand side of 
Eq. (|27|) . The function V ai ...a m in Eq. ( [27D is the m-parton dispersive contribution to the 



two-point function of the parent parton a. Being a two-point function, it is proportional to 
the unity matrix in the colour indices of the parton a. Thus, we can sum over the colours 
of the partons in the tree-level amplitudes, and we can rewrite Eq. (|2~7D in terms of the 
spin-polarization tensor 7^ introduced in Eq. (0): 

/o 2e \ m— 1 

M„ „.. (/>,. • • • ,p m , . . .)| 2 ~ t'bi + .-.+P™, ■ ■ •) K ( rL>i; -,Pmj 



Si. 



(29) 

The function Vj^ am is simply obtained from am in Eq. ( |2TD by performing the average 
over the colours of the parent parton a and extracting the factor in the round bracket 
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Figure 1: Dominant diagrams for the multiple collinear limit at 0(a™ ) in a physical 
gauge. 
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1 
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Figure 2: An interference diagram for the triple collinear limit a — > a±a2as. 



on the right-hand side of Eq. (|29|). Thus the tree- level function V^} (pi, . . . ,p m ) does 
not contain any additional power of the QCD coupling as- Note also that the spin tensor 
Xi"}(Pi + • • • + Pmi ■ ■ •) is n °t yet exactly the physical polarization tensor of Eq. ([25]) . In 
fact, the momentum of the parton a is off-shell ((pi + . . . + p m ) 2 = Si... m ^ 0) and, thus, 
Xi n ,(Pi + • • • + Pm, ■ ■ •) is gauge-dependent. 

To proceed, we should consider separately the two cases in which the parton a is either 
a quark (or antiquark) or a gluon. 



Quark splitting processes 

It is convenient to include the spin-polarization matrices p\ + . . . + j> m of the decaying 
quark a = q in the definition of the Dirac matrix Vj™ Qm (p\ , . . . ,p m ). The most general 
decomposition of V^} . 0m (pi, . . . ,p m ) is 

Vail.a m {Pii ■ ■ ■ iPm) ~ X/ (scalar amplitude) • (string of gamma matrices) . (30) 

Any string of gamma matrices is obtained by multiplying an arbitrary number of terms fa 
with I = 1, . . . , m + 1, where f can be either fa = j> i or ^ m+ i = iis\_ m jn ■ (pi + . . . + p TO ). 
The matrices ^, like V^ >0 (pi, . . . ,p m )> are homogeneous functions of n M with vanishing 
homogeneity degree. Thus, by Lorentz covariance, the amplitudes on the right-hand side 
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of Eq. (|30"D are scalar functions of the sub-energies s^- and the longitudinal-momentum 
fractions Zi = n-pi/n- (pi + . • -+p m )- Moreover, they are rational functions of the variables 
Sij, Zi and thus, by dimensional analysis, the corresponding strings can contain only an odd 
number of gamma matrices. 

We can now exploit the hermiticity properties of V^} Mm (pi, . . . ,p m ). Since the scalar 
amplitudes are real, the strings of gamma matrices appear in the form 

, l x(fc-l)/2 

) [1)ix1>n'--ti k +1>i k ---1>idh} > (31) 

where the normalization by the overall power of l/si... m has been introduced to make the 
scalar amplitudes on the right-hand side of Eq. (HO) dimensionless. Owing to the fact that 



k is odd, the terms with k = 3, 7, 11, . . . in Eq. ([H]) can in turn be reduced to strings 
that contain k — 1, 5, 9, . . . gamma matrices by using the anticommuting properties of the 
Dirac algebra. This is the simplest form in which we can write the general decomposition 
of Eq. 



We can now discuss separately the cases that involve the collinear decay of less or more 
than four partons. 

From the previous discussion we conclude that, when m < 3, the functions V^} Mm can 
be written as follows 

VgUtei, ■■■,Pm) = t 4 9 \{sfl, A + B^({s jh Zj }) - Jl l - n \ n ; , (m < 3) , 

i = l " ' \Pl ~r • • • ~r Pm) 



(32) 



while, when m = 4, we have 



4 



K ( rL 4 (pi,---,P4) = E^ g) (W^»A + g(g) (W^» { is '-\ (33) 

j =1 n ■ (Pi -r ■ ■ ■ ~r Pa) 

+ ° {{S ^ ZA) s 1 .. A n.( Pl + ...+p i )- (M) 

Then we can proceed to single out the dominant singular behaviour of Eqs. (|32p and ([53) 
in the multiple collinear limit. Since the scalar functions Af'({sji, Zj}), B^ q \{sji, Zj}) and 
G {{ s jh z j}) are dimensionless, they are invariant under the scale transformation fl23|) . 
Moreover, since we can write 

rf = z i (p 1 + ...+p m y + k? + 0{k? ± ) , (35) 
by rescaling the transverse momenta as in Eq. (^) we obtain the following scaling behaviour 

m 

V^.. a JPi,---,Pm) = (A + ---+An) E^4 (9) ({^,^})[1 + 0(A)] , (m<4). (36) 

i=i 

Thus, inserting Eq. ( ^B] ) into Eq. (E^), we can use the spin polarization factor j> x + . . . 
to reconstruct the matrix element squared |A^^.(Pi + ••• +Pm,---)\ 2 - Having already 
factorized the dominant singular term, we can now replace p\ + . . . + p m xp in | 2 , 
so that its gauge dependence disappears, and we obtain the factorization formula 
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Moreover, we also obtain an explicit expression for the quark splitting function in terms of 
the scalar amplitudes A^f 1 in Eqs 



and (B3): 



^E*4 tf) ({«ii>*i})> (m<4). 



a\ ...a r , 



(37) 



This argument to prove collinear factorization is based on the fact that a single spin 
structure (see Eq. (|36|) ) dominates the collinear limit of the quark decay function am . 
In particular, this implies that spin correlations are absent from the collinear decay of a 
fermion, independently of the number of its spin polarizations. However, the argument 
works only for the cases with m < 4. When m > 4 collinear factorization still applies 
but, as shown below, spin correlations cancel only if we use a dimensional-regularization 
prescription in which the massless fermion has two spin polarizations. 

According to our definition, the scalar amplitudes on the right-hand side of Eq. (|30|) 
are dimensionless and, hence, they are invariant under the scale transformation (^). The 
collinear limit of Eq. (pOf ) is thus determined by that of the strings of gamma matrices in 
Eq. (|3lD. Using Eq. d35|) and rescaling the transverse momenta as in Eq. (p3|) , we see that 
the strings that dominate in the multiple collinear limit are those of the form 



S1...T 



SI...T 



K + • • ■ An) &i ' ' ' ty'hk + ^*2fc " "$n fih \Pl + • • • Pm)) 

X \jf> fa fti 2 ■ • • fti 2k + ljti 2k " 'jti 2 fii-L fl\ [1 + 0{\)\ , 



(3? 



where the dots stand for the product of % factors. We can now multiply Eq. (|38|) by unity 



in the form 1 



+ njij)) j (2p ■ n), and, using {fa ^} = and f = 0, we obtain 



si. 



Si. 



X'jf) 



+ Tfip plfb + 



l2k 2p ■ n 

l f i 2 - " - f-i 2k i 1 + i 1 f-i 2 k 



2p-n ' rl2k 



2xp ■ n 



fti2k " ' fl2 fill V 

xp 1 . 



(39) 



Denoting by Xs{p) the spinor of an on-shell fermion with momentum p and spin polarization 
s, we then replace the polarization matrices xp 1 in Eq. (j39D by using the identity xp 1 = 
J2sXs(xp)Xs( x P) an d we can rewrite the string in Eq. ( j39D as follows 



[ Xs(xp) Xs'( X P) 



si. 



— ( \ fth fci2 ' ' ' l^»2fc j 1 j 1 &hk fife f-h 

Xs\ X P) o 

2xp ■ n 



Xs'(xp) . (40) 



When inserted in Eqs. (|30| ) and (p9|) , the factor in the square bracket reconstructs the po- 
larization matrix of the decaying quark and, thus, the spin-polarization tensor T q ss (xp, . . .) 
in the factorization formula (j25D . The remaining factor in Eq. fl40|) gives the contribution 
of the string in Eq. (|38D to the quark splitting function P^... am - 



By explicit construction we see that in general the splitting function p£ s ' am is not 
diagonal with respect to the spin indices. Nonetheless, the spin correlations are absent 
within the dimensional-regularization prescription used throughout the paper. Since we are 



10 



considering only two helicity states for massless quarks, we have Xs=±i(p) — |(1 ±75)x(p)> 
where x{p) is a generic Dirac spinor. Thus, using the general properties of the Dirac algebra, 
the contribution of Eq. (^) to the splitting function can straightforwardly be recast in a 
form that explicitly shows the cancellation of the spin correlations: 



5 s * 



k It- V? 



1 



Sl„, 



Tr 



2xp ■ n 

iji (g fj^ fkj 2 ■ ■ ■ f,j 2k + f,j 2k ■ ■ ■ f,j 2 fj^ g) 
Ap ■ n 



s, s 



±1 



(41) 



where Tr denotes the trace of the Dirac matrices. The identity in Eq. ([IT]) relies on the 
definition and the properties of the chiral projectors |(1 ±75) and the absence of spin cor- 
relations ultimately follows from helicity conservation in the quark-gluon vector coupling. 

This method to derive collinear factorization also provides us with an expression of the 
(spin-averaged) quark splitting function in terms of the dispersive part V^} am (pi, . . . ,p m ) 
of the two-point quark amplitude. From Eqs. (|30|) and (f4~UD, we find 



(-Pai...a m ) — Tr 



>K ( " } .a m (Pl>---^ 

4 n ■ (pi + . . . + 

Pm 



(42) 



This equation is useful for a straightforward evaluation of the splitting function for the 
multiple collinear decay of a quark. 



Gluon splitting processes 

Unlike the quark case, it is convenient to define the gluon two-point function V^ n \ m without 
including in it the spin-polarization tensors d^)(pi + ■ ■ ■ + p m ) of the two external gluons. 
Because of Lorentz covariance and the vanishing degree of homogeneity with respect to n M , 
the spin tensor V^f)a can be decomposed in terms of dimensionless scalar amplitudes as 

VaZ^iPU ■ ■ ■ ,Pm) = A^({8 jh Zj}) <T + £ B%({ Skh Z k }) ^ 

i,j=l s l...m 

+ V r {9) ( Is , z] ) P ^ + ^ I D [q] (\c , z \) — Sl - m (A3) 
Then, we have to multiply V^^} m by ^he g mon polarization tensors as follows 

€ H (P1 + ...+Pm) V^liPl, ■ ■ ■ ,Pm) d£ (n) (pi + ...+ Pm ) ■ (44) 

Inserting Eq. (|4"3| ) into Eq. (|44T), we immediately see that the scalar amplitudes C\ ({sji, Zj}) 
and D^({sji, Zj}) give a vanishing contribution because the gauge vector is orthogonal 
to the polarization tensors. As for the second term on the right-hand side of Eq. (pE5D, 
we can extract its dominant collinear contribution by simply performing the replacement 
pf — > fcf . Indeed, using Eq. (|35|) and 

(pi + . . . +Pm) fJ ,d^(p 1 + . . . +p m ) = C(si... m ) , (45) 
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we have 

€(n)(Pi + ---+Pm) tfi p] d; {n) ( Pl +...+ Pm ) = k? K + 0(X 3 ) , (46) 
so that the longitudinal component of is suppressed in the multiple collinear limit A — > 0. 

We can now safely perform the collinear limit and we obtain the factorization formula 
(531) and an explicit expression of the gluon splitting function in terms of the scalar ampli- 



tudes in Eq. fl43|) 



m k^k v 
PZ..a m = ^ 9 \{ Sjl , z 3 }) gT + £ B${{8 U , z k }) -±-L . (47) 

i,j=l s l...m 

The splitting function can be averaged over the spin polarizations of the parent gluon 
according to Eq. (f|), and we obtain 

1 1 171 , ^ h- ■ h ■ 



(P ai ...aJ = ^T^y dM PZ..a m = -A^{{S JU Zj})- £ B^({s kl , Z k }) 



Si... 



(48) 

where 

m 

2fcj • kj Sij ^ (ZiSjk + Zj&ik) + 2z$2jSi...j7i • (49) 
fc=l 

Note that, since 2ki ■ kj = 2 Pi ■ d( n )(pi + ■ ■ ■ +p m ) • ^(n)G°i + • • • +Pm) • Pj, the spin-averaged 
splitting function can also be expressed in terms of the Lorentz trace of Eq. 



(P ai ...aJ = ~ 2(1 l _ e) <(„)(Pl + ■ ■ ■ + Pm) V a \ P }2l( Pl , ..., Pm ) d pfl(n) ( Pl + ...+p m ) . (50) 

In the following subsection we present the explicit calculation of the quark and gluon 
splitting functions in the triple collinear limit. 



2.4 Collinear splitting functions at 0(a\) 

The list of (non-vanishing) splitting processes that we have to consider is as follows: 

s ->• si + S 2 + Ss , (s -> si + S 2 + Sa) , (51) 

S -> Si + S2 + S3 , (S -> Si + S2 + S3) , (52) 

S -»• Pi + P2 + S3 , (S -»• Pi + 92 + S3) , (53) 

9 Pi + S2 + S3 , (54) 

9 ^ 9i + 92 + 93 • (55) 

The superscripts in g', if denote fermions with different flavour with respect to q, q. As 
already mentioned in Sect. |2.3| , the splitting functions for the processes in parenthesis in 
Eqs. ( |51p and (|52|) can be simply obtained by charge-conjugation invariance, i.e. Pq' iq ' 2 q 3 = 

P q '^q 3 and Pq iq2 q 3 = P qi q 2 qz- m summary, we have to compute five independent splitting 
functions. 
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Figure 3: The diagram for the collinear decay q — ► q' x ^ 2 q^. 



To illustrate our calculation, we first consider the process in Eq. (0), that is, the 
case in which a quark-ant iquark pair <^g 2 and a quark q$ with different flavour become 
collinear. This is the simplest case, because the two-point function V^, q3 {pi,p2,Pz) for 
the corresponding splitting process is obtained by squaring the sole Feynman diagram in 
Fig. |||. According to the definition in Eqs. (p9| ) and (^), we extract the overall factor 
(87r/i 2e as/si23) and, performing the average over the colours of the decaying quark q, we 
find 



V£2 aga (pi,P2,P3) 



+ 



1 n rp S 123 

2 S12 



Zl + z 2 
2^12,3 



Zi + Z 2 )Si2 




Zl + Z 2 \S123 



j -S123 
(Pi +P2 +P3y 



(56) 



where t 12 ,3 is the kinematical variable defined in Eq. (22). Note that Eq. ([56]) has the 
general structure obtained in Eq. fl32"l) . Using Eq. (|42]) to compute the splitting function 
(Pq' iq ' 2 q 3 ), the last two terms on the right-hand side of Eq. fl5"6]) give a vanishing contribution 
and we obtain the final result: 



^ r< rp S 123 

— O pi ;? 

2 F * S12 



*12,3 + 4 ^3 + (*1 



^2 



s 12 s 123 



*L + z 2 



+ (l-2e) (zi + z 2 - 



S12 
S123 



(57) 



The calculation of the splitting functions for the other processes in Eqs. (|52])-(p5]) can 
be performed exactly in the same manner, by using the general procedure discussed in 
Sect. pT3"[ We first compute the corresponding two-point functions V^} 2a3 (pi,p 2 ,p 3 ) and 
then, using Eqs. (pE2"|) and P7j), we evaluate the splitting functions. Since the intermediate 
expressions for the two-point functions are quite cumbersome, in the following we limit 
ourselves to showing the relevant Feynman diagrams and to presenting the final results for 
the splitting functions. 

The calculation of the splitting function for the case of final-state fermions with identical 
flavour involves a diagram analogous to that in Fig. ||]plus its crossed diagram (see Fig. |j). 
Thus, the result can be written in terms of that in Eq. (|57j), as follows 



9l<32<?3/ 



3) + (Pqiqlq 3 ) 



(51 



13 



Figure 4: The diagrams for the collinear decay q — > q^qz- 



where the interference contribution is given by 



3(id) \ 
919253/ 



+ 





'c F - 


-Hf' 1 




Sl23 


"1 + 


zj 2z 2 




Sl2 


1 - 


Z2 1-23 




S 123 Z l 


[ 1 + *? 




susn 2 


(i-*0(i- 


~*0 



/2S23 _ 



1 — Z2 1 



2z 2 



*3. 



! (l-^3) 



el 



1-^3 



+ (2^3) . (59) 



{ 



9 2 



93 

* — — > + 





(1<->2) } 




Figure 5: The diagrams for the collinear decay q — > 5 , i5 , 2<?3- 



The splitting function of the remaining quark-decay subprocess is obtained by squaring 
the diagrams in Fig. ^. It can be decomposed according to the different colour coefficients: 



( P 9l92qs) - C F ( P gig 2 q 3 ) + C F C A (Pg\gtq 3 ) 



(60) 



and the abelian and non-abelian contributions are 



/ p(ab) \ 
\ .91.9293/ 



+ 
+ 



'123 



2si 3 s 23 



^3 



1+j 

^2 



e^±^-e(l + e) 



2^2 



Sl23 



z 3 (l-zi) + (I-22) 



Sl3 

(1-0 



^1^2 



+ e 2 (l + 23) - e(^ 2 + z x z 2 + z 2 



1-^2 
^1^2 



( £23 

S13 



+ (1 - 2) 



(61) 
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/ p(nab) \ 
\ 919293' 



+ 
+ 
+ 



.1- €)[-?£■ + ---) + 



2 

123 



4s( 2 '4 2) 2s 12 s 13 



;i-z 3 ) 2 (l-e) + 2z 3 



z 2 2 (l-e) + 2(l-z 2 ) 



l-z 3 



'123 



4513^23 



23 



Z 2 

[l-z,ni-e) 



2z, 



Z\Z 2 



+ 6(1-6) 



g 123 
2^2 
■S123 

2s 13 



^ _ c ^i(2 - 2gi + z|) - z 2 (6 - 6z 2 + z|) 2e ^3(^i - 2z 2 ) - z 2 



1-e) 



^2(1 - 23) 

'3 



;i-z 2 )3 + zf _z 2 /2(l-Z 2 )(z 2 -Z 3 ) 



^2(1 - z 3 ) 
:,(! - zQ + (1 - z 2 ) 3 

ZlZ 2 



22(1 - 23) 



22(1 - 23) 

- Zl + z 2 \ 



z? + z| 



ZlZ 2 



+ (!<-> 2) . 



(62) 



As discussed in Sect. p73j in the case of collinear decays of a gluon (see Eqs. (fjj, pq)), 
spin correlations are highly non-trivial. 





Figure 6: The diagrams for the collinear decay g — > gi<?2<?3- 

To compute the splitting function for the decay into a qq pair plus a gluon, we have 
to evaluate the square of the diagrams in Fig. ||. The colour-factor decomposition of the 
splitting function is 



pfiu _ p rp fin" (ab) . p rp A/if (nab) 
9i<?2?3 ^F-LR r g 1 q 2 q-i "r /? 919293 

where the abelian and non-abelian terms are given by 

2Sl23 s 23 + (1 ~ e )( s 123 ~ s 23) 2 



P; 



/zf (ab) 



919293 



+ 



-5 

43123 
S12S13 



-2 + 



^A^^ ~\~ k- 



S12S13 

1 - e)kfk» 



(63) 



(64) 
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919293 
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Sl23 



S12S13 
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S123 



z 3 
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S12 Zi(l - Zi) S 23 Zi(l - Zi) 
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+ 



gl23 
S12S23 



-2s 123 g 



lxv z 2 {l-2z 1 ) 

*i(l 
( 2z 2 (z 3 



\ Zi(l - Zi) 



X6 k^ 



+ (1 




3) 



(65) 




Figure 7: The diagrams for the collinear decay g — ► g\g 2 gz- 



In the case of gluon decay into three collinear gluons we have to consider the diagrams in 
Fig. [7]. Note that the contribution of the four-gluon vertex cannot be neglected. The result 
for the splitting function is quite involved. Its expression can be simplified by exploiting 
the complete symmetry under the six permutations of the gluon momenta. We obtain 
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+ (5 permutations) . 



(66) 



The splitting functions in Eqs. 
the parent gluon according to Eq. ( 

Performing the average, we obtain 

,p(ab)_ v 
\ r giq2q3/ 



2(1 -c) 



can be averaged over the spin polarizations of 



2 - (1 - e)s 23 ( — + —)+ 2 6 
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(69) 



V 919293/ 



+ 



+ 




2i^2 — 1 2i^2 — 2 



1 - 2:3 ^3 

2 2 )(l- 22 3 ) 



2 + 2 Z3 



2 3 2i(l - 2 X ) 
l+2^i(l + ^i) 
2(l-2 2 )(l-2 3 ) 



^3(1 - Zs) 



+ 2 2 2 3 - 2 + 



z 1 (i + 2z 1 ; 



1-221(1-20 

22 2 2 3 



+ (5 permutations) 



(70) 



The (9(ag)-collinear limit of tree- level QCD amplitudes has been independently consid- 



ered by Campbell and Glover | 2T| . They have computed only the spin-averaged splitting 
functions. The comparison with their results has been discussed in detail in Ref. |22] and 
we do not repeat it here. Our results agree with those of Ref. [^TJ. Since the methods used 
by the two groups are completely different (cf. the discussion in Sect. [I]), this agreement 
can be regarded as an important cross-check of the calculations. 

The expressions of the spin-dependent splitting functions Paf a2 a 3 derived in this section 
refer to the CDR scheme. Other dimensional-regularization schemes can be used. We 
mention two of them, which differ from CDR only by the number of spin-polarizations of 
quarks and gluons. 



The dimensional-reduction (DR) scheme |31[ works by considering two spin-polarization 
states for quarks and two for gluons. The corresponding spin-dependent splitting functions 
P^lcaaz are srm ply obtained from those in the CDR scheme by setting e = 0. 

The 'toy' dimensional-regularization (TDR) scheme introduced in Ref. [[Kj considers 
d— 2 = 2(1— e) spin-polarization states for quarks as for gluons. Its practical implementation 
is very simple. When computing traces of gamma matrices, we should use the relation 
Tr 1 = 4(1 — e), where 1 is the unity matrix in the spinor space. The corresponding spin- 
dependent splitting functions P^ a2a3 are obtained from those in the CDR scheme by the 
simple replacement T R — > T R (1 — e). 

The QCD results presented in this section can also be extended in a straightforward 
way to the abelian and supersymmetric cases. 

In the case of QED, we have to perform the replacement as — > ot in the factorization 
formula (|25|), and the relevant splitting functions, P^^a 3 \ f° r the triple collinear limit are 
obtained from the QCD splitting functions as 



piQED) 



2 2 



,(P, 
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p(QED) Up \ {n) 

1 519293 q\ 919293 y ab ' \ ts -l 

p(QED) = 4/p 

717293 9 V 9 

p(QED) A(p 

71 9293 9V ' 



919293 
919293 



ab. 
ab. 



where e g is the quark electric charge and the notation (• • •) b means that we have to set 
Cp = Tr : = 1 and Ca = in the QCD expression inside the round bracket. 

The supersymmetric version of QCD, namely N — 1 supersymmetric Yang-Mills theory, 
is obtained by replacing the quark with the gluino g, a Majorana fermion in the adjoint 
representation of the gauge group. To obtain the corresponding splitting functions, we have 
to change the colour factors accordingly, and we have to identify q = q = g after having 
summed over the different permutations of the final-state fermions. We have 

P~ - - = (P- 4- P - + P - ^ 

1 g\gigz y- 919293 ~ 1 919293 ~ 1 919293 /sqcd ' 

^919293 = (^19293) g QCD , ( 72 ) 

-^919293 = (-^919293 -^919293) SQCD ' 

where the notation (...) means that we have to set Cp = 2Tr = Ca in the QCD 

V / SQCD 

expression inside the round bracket. 

Gluino and gluon amplitudes are related by supersymmetry transformations. In the 
collinear limit, these transformations relate the total splitting functions Pjjl+3 and Pjj^s for 
gluino and gluon decays, which are defined as 



pss _ pss _|_ 

g~*3 — ^919293 ' 



PgLss + (3-1) + (3-2) = S ss ' (P^ 3 ) , (73) 



= P 9 r 9293 + K293 + (1 ~ 2) + (1 <- 3)j . (74) 
In the four-dimensional supersymmetric theory, gluon and gluino have the same decay 
probability. Provided supersymmetry is not broken by the dimensional-regularization pro- 
cedure, we thus have the following supersymmetric Ward identity: 

P^ 3 = -ST <^3> • (75) 

Note that the Ward identity holds for the spin-dependent splitting functions. Since spin 
correlations are absent in the gluino splitting function, they cancel in the right-hand side 
of Eq. fl75D ) and P|^3 an d Pg^>3 differ only by the overall spin-factors 5 SS and —g^ u . As is 
well known, the splitting functions PJ^ an d are related by a similar Ward identity 
at 0(a s ). 

The identity (|75|) is violated in the CDR scheme, because gluinos and gluons have a 
different number of spin-polarization states. The Ward identity is recovered in the e — > 
limit or, equivalently, in the DR scheme, which is known to explicitly preserve supersym- 
metry. Our results for the spin-dependent splitting functions fulfil Eq. ([75]), and this is an 
important check of our calculation. As pointed out in Ref. [Q, the Ward identity at 0(as) 
is fulfilled also in the TDR scheme. We have verified that this remains true at C?(a|), as 
expected from the fact that in the TDR scheme the number of gluino states is the same as 
the number of gluon states. 
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3 The soft behaviour 



The tree-level matrix elements M.(px,P2, ■ ■ • ) are singular not only when parton momenta 
become collinear but also when one or more of them become soft. In QCD calculations of 
physical cross sections at NLO, the soft limit is approached when the momentum of a single 
gluon vanishes. At NNLO we have to consider three different types of soft configurations: 

• the emission of a soft quark-antiquark pair, 

• the emission of two soft gluons, 

• the emission of a soft gluon and a pair of collinear partons. 

The behaviour of the tree-level matrix elements in these singular limits is considered in this 
section. We also discuss the generalization of the corresponding factorization formulae to 
higher perturbative orders. 

3.1 Colour correlations and eikonal current at 0(as) 

The emission of a soft gluon does not affect the kinematics (momenta and spins) of the 
radiating partons. However, it does affect their colour because the gluon always carries 
away some colour charge, no matter how soft it is. Unlike the case of soft-photon emission 
in QED, soft-gluon emission thus does not factorize exactly and leads to colour correlations. 

To take into account the colour structure (as well as the spin and flavour structures), it 
is useful to introduce a basis {\ci, c n ) ® |s 1; s n )} in colour + helicity space in such a 
way that the tree- level matrix element in Eq. (p]) with n final- state partons can be written 

as 



Ml\'-^ sl '---' Sn (p 1 ,...,p n ) = ((c 1 ,...,c n \®(s 1 ,...,s n \) \M ai ,..., an (p 1 ,...,p n )) . (76) 



Thus \M. ai ,...,a n (pi, ■ ■ ■ iPn)) is a vector in colour + helicity space. 

According to this notation, the matrix element squared (summed over final-state colours 
and spins) \M\ 2 can be written as 



To describe the colour correlations produced by soft-gluon emission, it is convenient to 
associate a colour charge Tj with the emission of a gluon from each parton i. If the emitted 
gluon has colour index c (c = 1, N 2 — 1), the colour-charge operator is: 



\M au .„ tan (p 1 ,...,p n )\ 2 = (M ai ,...,a n (Pl,---,Pn)\M 



(Pl,...,p n )) • 



(77) 



T t = (c\ If 




and its action onto the colour space is defined by 



Clj . . . , Cj, . . . , C m , C| J. j| Oi, . . . , Oj, . . . , 



b m ) — 5 Clbl ...T^ bi ...5 Cmb . 



(79) 
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where T£ b = if ca b (colour-charge matrix in the adjoint representation) if the emitting parti- 
cle i is a gluon and T% = t^g (colour- charge matrix in the fundamental representation with 
a, P — 1, ... , N c ) if the emitting particle i is a quark (in the case of an emitting antiquark 

rpa — Ja j.a \ 

1 a/3 — L a/3 ~ l f3a)- 

The colour-charge algebra isfj]: 

7fT? = Ti ■ Tj = Tj ■ Ti if i ± j; T\ = Q, (80) 

where Cj is the Casimir operator, that is, Cj = Ca = N c if i is a gluon and C{ = Cf = 
(N% — 1)/2N C if i is a quark or antiquark. 

Note that, by definition, each vector \M.ax,...,a n {Pu ■ ■ ■ ,Pn)) is a colour- singlet state. 
Therefore colour conservation is simply 

J2T l \M ai _ an ( Pl ,..., Pn ))=0. (81) 



Let us now consider the tree-level matrix element M. g , ai ,...,a n ( ( l,Pi, ■ ■ ■ ,Pn) i n the limit 
where the momentum q of the gluon becomes soft. Denoting by c and \x the colour and spin 
indices of the soft gluon, the matrix element fulfils the following factorization formula || 

(c; fl\ Mg (q, Pl ,...,p n )) ~g s ^J c ^{q) \M 

ai,...,a n ( P l,...,Pn)) , (82) 

where | M. ai ,...,a n {Pii ■ ■ ■ ,Pn)) is obtained from the original matrix by simply removing the 
soft gluon q. The factor J^(q) is the eikonal current 

J^q)=^T t -^-, (83) 

which depends on the momenta and colour charges of the hard partons in the matrix 
element on the right-hand side of Eq. ([82]). The symbol £ ~' means that on the right-hand 
side we have neglected contributions that are less singular than 1/q in the soft limit q —>■ 0. 
Note that Eq. (|82]) is valid in any number d = 4 — 2e of space-time dimensions, and the 
sole dependence on d is in the overall factor fi e . 



The factorization formula (82) can be derived in a simple way by working in a phys- 
ical gauge and using the following soft- gluon insertion rules. The coupling of the gluon 
to any internal (i.e. highly off-shell) parton in the amplitude M. g> ai > ...,a n { ( l-iPu ■ ■ ■ ,Pn) is 
not singular in the soft limit; it can thus be neglected. The soft-gluon coupling to any 
external or, in general, nearly on-shell parton with colour charge T and momentum p can 
be factorized by extracting the contribution gsfi e 2p' 1 T for the vertex and the contribution 
I/O 9 + — 1/G° 2 + 2p • q) for the propagator. 

An important property of the eikonal current is current conservation. Multiplying 
Eq. ( p2|) by we obtain 

q^(q)=f2Ti, (84) 
i=i 

"More details on the colour algebra and useful colour-matrix relations can be found in Appendix A of 
Ref. §. 
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and thus 



q^J fl (q)\M ai ,..., an ( Pl ,...,p n )) = £T 4 \M ai ,...,a n {pi,...,Pn)) = , (85) 



i=l 



where the last equality follows from colour conservation as in Eq. (|8ip . 



Although Eq. ( [82] ) is most easily derived in a physical gauge, the conservation of the 
eikonal current implies the gauge invariance of the squared amplitude summed over the soft- 
gluon polarizations. Squaring the eikonal current and introducing the gluon polarization 
tensor d^ u (q) = {—g^ u + gauge terms) in Eq. ([|), we have 

[J^rf cW?) J v {q) = - £ Ti ■ Tj ^ - q) + . . . , (86) 

where we have used the fact that the gauge terms in d fJiU (q) are due to longitudinal polar- 
izations proportional either to q^ or to q u . Thus, the dots on the right-hand side stand 
for gauge-dependent contributions that are proportional to the total colour charge J27=i Ti 
and, hence, that cancel when they are inserted in | Ai ai ,...,a n (pi, ■ ■ ■ ,Pn))- 

Using Eq. (|86|), the soft-gluon factorization formula at 0(azs) for the squared ampli- 
tude is 

n 

\M g , ai ,... ian (q, P l, ..., Pn )\ 2 ^ -47T« S /i 2e £ <%(?) \M^fl >an ( Pl , . . . , Pn )\ 2 , (87) 

i,J=l 

where the scalar eikonal function Sij(q) for the emission of a single gluon can be written in 
terms of two-particle sub-energies Sy = (pi + Pj) 2 as follows 



<%(g) 



Pi ' Pj _ 2Sjj 



(Pi ■ q) (Pj ■ q) Siq Sjq 



The colour correlations produced by soft-gluon emission are taken into account by the 
square of the colour-correlated tree-amplitude |A^^| 2 on the right-hand side. This is 
defined by 

\M^l an ( Pl , ...,Pn)\ 2 = (M ai ,..., an ( Pl , . . . ,p n ) | Ti ■ Tj I Ma u ... >an (pi, ■ ■ ■ ,Pn)) (89) 



M^^ipu . . . , Pn )\ v;:;,, •/;;;„, M^ : ^{ Pl , ..., Pn ), 

where the sum over the spin indices is understood. 

The soft-gluon factorization formula is often presented |2l]] in an equivalent way by 
decomposing the matrix element in terms of colour subamplitudes In this formalism, 
the eikonal function Sij(q) in Eq. (|88|) controls the factorization properties of the square of 
the colour-connected subamplitudes. 



3.2 Emission of a soft qq-pmr 

We now consider the tree-level matrix element M q ,q, ai ,...,a n (qi, <?2,Pi, • • • ,Pn) when the mo- 
menta qi and q2 of the quark q and the antiquark q become soft (q\, q<i — > at fixed q\jq-i). 
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In this limit the matrix element squared has the dominant behaviour: 



\M q> q >ai ,..., an (qi,q2,Pl, ■■■,Pn)\ 2 ~ 7 w 7 j 777 7 j 77 • (90) 

{Qi ■ Q2) [Pi ■ (9i + 92)] [Pj ■ (9i + 92)] 

When integrated over the phase space of the quark- ant iquark pair, this behaviour gives rise 
to a single-logarithmic soft singularity, in addition to possible single- and double-logarithmic 
collinear singularities. 

The soft singularity arises when the gg-pair is produced by the decay of a gluon that 
carries the soft momentum q\ + q 2 (Fig. ||). Thus, using the soft-gluon insertion rules 
described in Sect. |3.1|, we can straightforwardly derive the following factorization formula: 



\M q) q >ait ... jan (qi,q2,Pi, ■ ■ ■ ,p n )\ 2 - (4vr/i 2e a s ) 2 

• (M ai ,...,a n (pi,---,Pn)\ I(qq)(qi,q2) \M ai ,...,a n (pi,---,Pn)) ,(91) 

where 

I(ri)(Qi, 9 2 ) = [J\qi + q2)fd x ^ qi + q 2 ) W v ( qi , q 2 ) d vp {q x + q 2 ) J p {q x + q 2 ) (92) 

= [^(91 + 9 2 )] t n^(gi, q 2 ) J v {q^ + <fe) + • • • • (93) 

The insertion operator I( qq )(qi, 92) depends on the colour charges of the fast partons 
ai, . . .,a n and it is given in terms of the soft-gluon current J^(qi + 92) in Eq. ( |83"1) and 
of Il Atl/ (g 1 , q 2 ), which is the gg-contribution to the discontinuity of the gluon propagator: 

n^(9i, 92) = — ^ { -gTqi ■ 92 + 9^ + 9 2 M 9i } • (94) 

The dots on the right-hand side of Eq. ( |9"3"D denote the gauge-dependent contribution 
to the insertion operator I( gq -)(qi,q2)- This term is due to the longitudinal polarizations 
(proportional to (91 + 92)° or (gi + q 2 ) 13 ) of the polarization tensors d a p(qi + 92) in Eq. Q9~2|). 
Since Il >J,u (qi,q 2 ) is transverse (gi^IP^gi, q 2 ) = 0) and the soft current J^(qi + q 2 ) is con- 



served (see Eq. (84)), the contribution of the longitudinal polarizations is either vanishing 
or proportional to the total colour charge of the fast partons. Because of colour conser- 
vation (see Eq. (|3lD), we thus conclude that the gauge- dependent part of I^(qi,q 2 ) does 
not contribute to the factorization formula (pTj) . 

Inserting Eq. (^) into Eq. ([H]) and performing the Lorentz algebra, we obtain the final 
factorization formula 



\M q ,g, ai> ... } a n (q 1 ,q 2 ,p l ,...,p n )\ 2 ~ (47r/i 2e a s ) 2 T H £ Tij(qi,q 2 ) |A*£f.. >0n (pi, . . . ,p n )\ 2 , 



(95) 

where \M.%fl „ | 2 is the colour-correlated tree-amplitude of Eq. (|89|) and the soft function 
^(91,92) is given by 

T / \ (gt • 9i) (Pj ■ 92) + {pj ■ 9i) {Pi ■ 92) ~ (gj ■ Pj) (9i • 92) (Q(] , 

vMl)q2) (9i • 92) 2 \pi ■ (9i + 92)] ]Pi ■ (9i + 92)] 1 J 

_ 2(pj ■ Pj) (91 • 92) + [p, • (9i ~ 92)] [Pj ■ (9i ~ 92)] + , g7 x 
2(9i • 92) 2 [Pi ■ (9i + 92)] [Pj ■ (9i + 92)] 
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Figure 8: Soft-gluon insertion diagram for the emission of a soft qq pair. 



Note that both the expressions (0) and ( j97D can equivalently be used to compute Eq. fl95|) . 
The difference between the two expressions, denoted by the dots on the right-hand side of 
Eq. (^) , gives a vanishing contribution to the factorization formula ([)5|) because of colour 
conservation (see Eq. ([Si])). 



3.3 Soft current for double gluon emission 



The limit of QCD tree-amplitudes when the momenta of two gluons simultaneously become 



soft was independently studied by Berends and Giele [pl| and by one of the authors [ 20] . 
The singular behaviour of the matrix elements can be described in terms of factorization 
formulae given in terms of process-independent two-gluon currents acting either on colour- 
ordered subamplitudes |B| or on the colour space of the hard partons |2(J . 



The formalism of the colour subamplitudes was used in Ref. |2TJ] to derive explicit soft- 
gluon factors for the square of colour-connected and colour-unconnected subamplitudes. 
In this section we recall the formalism and the results of Ref. p0[ and we present the 
corresponding factorization formula for the square of the matrix elements. 



We consider the tree-level matrix element M. 



.,a„(9i,<?2, ,Pi, ■■■,Pn) when the mo- 



menta q\ and q2 of the two gluons become soft. The limit is precisely defined by rescaling 
the gluon momenta by an overall factor A: 



qi -» Xqi , q 2 -> Ag 2 , 
and then performing the limit A — > 0. The matrix element thus behaves as 



M 



9,9>- 



0(1/ A 2 ) + 



(91 



(99) 



where the dots stand for less singular contributions as A 
evaluating the dominant singular term Oil/ A 2 ). 



0. We are interested in explicitly 
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Note that the double soft limit is more general (accurate) than the soft limit in the 
strong-ordering approximation, that is, when p$ ^> q\ ^> q 2 . The strongly-ordered limit 
describes only the double-logarithmic soft singularity of the matrix elements. The double 
soft limit reproduces consistently the double-logarithmic behaviour and correctly evaluates 
also the single-logarithmic soft singularity. 

We denote by ai,a 2 and //i,/U 2 the colour and Lorentz indices of the two gluons, re- 
spectively. In the double soft limit the matrix element fulfils the following factorization 
formula |20| 



(ai, a 2 ; fi X , n 2 \ M gi g iai ,..., an (qi, &,pi, ■ ■ ■ ,p n )) - 9 S V e J «/S(9i> <Zt) I M ai ,...,a n (pi, ■ ■ -,Pn)) , 

(100) 

where the two-gluon soft current J^liqi, 92) is the generalization of the eikonal current in 
Eq. (H). 



The explicit expression of the soft current is [20 



« 2 0?1^) 



„M1 „M2 
rpai fi rpa 2 rj _|_ 



+ E 



I* Pi - f /i ' ■ ( h 

a±a rpai ri 



S aia T? 



Pi ■ Q2 



M2 

^ ja2aia "1 | rpa 



Ml 

Pi 



Qi -Q2 J 1 Pi ■ (qi + Q2] 



+ 



_|_ I fia 2 a rpai Pi ^ja\aia ^2 



Ml 



rpa 



p 



M2 



Pi ■ Qi 



Qi ■ Q2J 1 Pi ■ (91 + 92) 



+ 



+ l i faaxa 2T a 9^ Pi \ (gg ~ gl) 

2 1 qi ■ q 2 Pi ■ {<h + Qi) 



(101) 



It can be derived by working in a physical gauge and using the soft-gluon insertion technique 
described in Sect. [Q[ We have to consider the diagrams in Fig. ^. The contribution on 
the first line of Eq. fllOlD comes from the eikonal emission of the two soft gluons from 
two different external partons (diagrams (a) in Fig. ^). The first term on the second and 
third lines come from the eikonal emission of the two gluons from the same external parton 
(diagrams (b) in Fig. |^). The remaining contributions in Eq. ( |101[ ) are proportional to 
faa 1 a 2 an d originate from the non-abelian diagrams of Fig. |9| (c). Note that the three-gluon 
vertex has to be treated exactly, without introducing any soft approximation. 




+ (l<->2) + 



)] 




(a) 



(b) 



(c) 



Figure 9: Soft-gluon insertion diagrams used to evaluate the two-gluon current J^laHqi, 92) ■ 
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The two-gluon current in Eq. ( |101| ) can be recast in the following equivalent form 



« 2 2 (9i,g 2 ) = \ {J£(qi) , ■/£(<&)} (102) 
pfgf -pf g^ 1 P%-{qi-q2) 



8=1 



(?i • ^2) [Pi • (ft + <?2)] 2[pi • (gi + g 2 )] 



(Pi ■ qi){pi ■ Q2) gi-92 



where the first term on the right-hand side is the colour anticommutator of the single-gluon 
eikonal currents of Eq. (|83D . This is the only contribution that survives in the abelian case, 
where it reduces itself to the product of two independent single-gluon currents. The second 
term on the right-hand side is typical of the non-abelian theory. 



Note that, as in the single-gluon case, the expressions ( |101| ) and ( |102| ) for the two-gluon 



current do not explicitly depend on the number d = 4 — 2e of space-time dimensions. 
However, because of the contribution proportional to g^ 1 ^ 2 in J^i 2 , an explicit dependence 
on the number d — 2 = 2(1 — e) of gluon polarizations appears (see Eqs. (|108|) and (|109|) ) 
by squaring the factorization formula ( |100| ). 

The current ^^2(91,92) fulfils the following properties. 

• It is symmetric under the exchange of the two soft gluons, 

^(9i,92) = ^(92,9i) • (103) 

• Its divergence is proportional to the total colour charge of the hard partons: 

gf -Oft. «0 = ^(92) s aia + -f aia2a -^- £i; a , (104) 

V 1 9i • 92/ 



i=i 



q^j;\Z{qiA2) = [j;\{qx)Ka + y a ^a J ^)Y, T t ■ (105) 

V 1 9i • 92/ i=1 

This property is analogous to Eq. (|4]) for the single-gluon emission and implies that 
the two-gluon current is conserved when it acts on a colour singlet state: 

9^^(91, 9 2 )| M au ...(p u ■■ •)) = 92 M2 ^(g 1 ,g 2 )|^ ai ,...(pi,...)) = 0. (106) 
Thus, the factorization formula ( |100| ) is gauge- invariant. 



In the strong-ordered limit q 2 <C gi, the third and fourth lines in Eq. (|101|) give 
subleading contributions and the current becomes 

« 2 (9i,9 2 ) - (jE(fc) 5 aia + if aia2a ^-\ jr(gi) • (107) 
\ 9i • 92 / 

Thus, the current correctly factorizes into the product of the two eikonal currents cor- 
responding to the iterative application of the leading-order factorization formula 



The double soft limit of \M. g ,g, au ...,a„(qi, 92, ,Pi, ■ ■ ■ ,Pn)\ 2 is obtained by squaring Eq. ( |100| ) 
and by summing over the soft-gluon polarizations. The square of the two-gluon current 
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involves a quite cumbersome colour algebra. Nonetheless, we find that the final result can 
be recast in a relatively simple form: 

J% a2 (qi, 9 2 )] f dr{q x ) d^(q 2 ) J% a2 (q 1} q 2 ) = \{j\q{) ,J 2 fe)} 

n 

- C A Y, T i- T i SiMii «0 + ■ ■ ■ , (108) 

where, as in Eq. ([86]), the dots stand for gauge-dependent terms. These are proportional to 
the total colour charge of the hard partons and, thus, give a vanishing contribution when 
inserted on | A4 ai ... 0n 

(p 1 ,...,p n )). 



The first term on the right-hand side of Eq. (|108|) is the only one that survives in the 
abelian case. It is given in terms of the colour anticommutator of the squares of the single- 
gluon currents in Eq. (|86|) . The second term is proportional to Ca and, hence, is purely 
non-abelian. It is given in terms of the two-gluon soft function Sij(qi, q 2 ): 



Sij(qi,q 2 ) 



[1 - e) pi - q x Pj ■ q 2 +Pi-q2 Pj • 9i 



(?i • Q2) 2 Pt ■ (?i + 92) Pj ■ (Qi + Q2) 



2 Pi ■ <?i Pj -Q2Pi- 92 Pj ■ q\ 



2 - 



Pi ■ qiPj ■ q2+Pi- q2Pj ■ qi 
Pi ■ (9i + 92) Pj ■ (91 + 92) 



(109) 



+ 



Pi ■ Pj 



+ 



29i ■ 92 [Pi ■ 9i Pj ■ 92 Pj -q\Pi- 92 



Pi ■ (9i + 92) Pj ■ (9i + 92) 
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(Pi-QiPj ■ Q2+Pi-q 2 Pj ■ qif 
Pi ■ 9i Pj -q 2 Pi- 92 Pj ■ 9i 



Expression (|109|) can also be written as 



c / v ( B .o.)/ n . Pi-Qi Pj ■ <l2 + Pi- q2Pj ■ 9i 

Oij{qi,q 2 ) = S.j 9i,92 + 7 — r j : r 

Pi ■ (9i + 92) Pj ■ (9i + 92) 

2pi ■ pj 



1_£) -Ut } M 



(9i ■ 92 



9i • 92 Pi ■ (9i + 92) Pj ■ (9i + 92) 



(110) 



where S^' '^ is the approximation of the soft function Sij(qi, q 2 ) in the strong-ordering limit 
(either q\ <C q 2 or q 2 <C qi): 

Pi-Pj ( 1 



(Pi ■ Pj f 



s$"->fa,<h) = 2-2- : + = 

9i"92 \Pi ■ 9i Pj ■ 92 Pj ■ 9i Pi ■ 92 / Pi- 9i Pj -QsPi' 92 Pj • 9i 



(111) 



Using Eq. ( |108j ), we can write the soft-gluon factorization formula for the square of the 
matrix element as follows: 



|-M 9 , 3 ,a 1 ,...,a„(9i,92,Pi,---,Pn)| 2 ^ (4vra s /i 2e ) (112) 

in n 

E <^(9l) 5„(fe) |^!i?(Pl. • • • ,Pn)| 2 - C A E <S^( gi , 92)|^£f. ) ., a J 2 



i,j,k,l=l 



i,i=l 
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where Sij(q) is the soft function in Eq. d88|) and a P ^ s ^ ne colour-correlated ampli- 

tude in Eq. (p9|). We can see that the double soft limit involves colour correlations that are 
more cumbersome than those appearing in the case of single-gluon emission. Indeed, the 
amplitude \M^} k ^\ 2 on the right-hand side of Eq. ( |112| ) is defined by 



l>f£?.!5?(Pi. • • • iPn)\ 2 = (M ai _ an (Pu ■ ■ ■ ,Pn)\ {Ti ■ Tj , T k ■ T,} \M au ..., an ( Pu ■ ■ .,Pn)) , 

(113) 

and leads to irreducible correlations among four different hard partons. 

The results discussed in this subsection can be presented in a different manner by using 
the colour subamplitude formalism. Considering the projection of Eq. ( |100| ) onto colour- 



ordered subamplitudes, it is straightforward to check that the colour current Jal^ijixi 92) 
leads to the colourless current derived by Berends and Giele (see Eqs. (3.11) and (3.18) 
in Ref. The square of this colourless current is denoted by S iqiq2 j in Sect. 5.3 of 



Ref. |2T| and is related to the soft function <%(gi, (72) in Eq. (P-09Q . More precisely, using 



the following relation 

n 1 ra 

T t ■ Tj [S i:j {qi, q 2 ) + S y ( 9l ) <%(g 2 )] = - £ T, ■ Tj S mq23 + . . . , (114) 

i,j=l i,j=l 

the second term on the right-hand side of Eqs. ( |108| ) and ( |1 1 2j ) can equivalently be written 
in terms of Si qiq2 j. The contribution denoted by the dots on the right-hand side of Eq. ( |114j ) 
is proportional to the total colour charge of the hard partons and, thus, it vanishes when 
inserted in the factorization formula ( |112| ). 



3.4 Soft— collinear limit at 0{a\) and at higher orders 

We now consider the tree-level matrix element Ai giai) ... jan (q, Pi, ■ ■ ■ ,p n ) in the limit where 
the momentum q of the gluon becomes soft (q — > 0) and, at the same time, two partons, 
say pi and become collinear. The collinear region is parametrized as in Eq. (B) and we 
are interested in the limit k±_ — > 0. 

Studying this soft-collinear limit we can neglect i) contributions that are uniformly of 
0(q) when q — > 0, and ii) contributions that are uniformly of 0(k±) when k± — > 0. The 
terms in class 1) are not singular in the soft limit and their contribution in the collinear 
limit can thus be taken into account by supplementing the results of this section with the 
(9(as)-collinear factorization discussed in Sect. |2.1| . Analogously, the terms in class ii) are 
not singular in the collinear limit and their contribution in the soft limit can be taken 
into account by supplementing the results of this section with the soft-gluon factorization 



formula at O(as) presented in Sect. IO 



This comment can be summarized in a formal manner by writing the square of the 
matrix element as 

\M g , ai ,a 2 ,...,a n ( ( l,Pl,P2,---,Pn)\ 2 = " F(q, p U p 2 , ■ . . , Pn) ■ (H5) 

Sl2S\ q S 2q 

The first factor on the right-hand side contains the correct scaling behaviour in the soft and 
collinear regions. Thus, the soft-collinear limit is defined by the soft {q — > 0) and collinear 
(k± — > 0) approximations of the function F(q,pi,p 2 , . . . ,p n ) at fixed ratio q/kj_. 
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To compute the soft-collinear limit we perform first soft approximations and then 
collinear approximations. 

The singular behaviour of the matrix element in the soft limit (and at fixed q/k 2 ^) is 
given by a factorization formula analogous to Eq. ( [52]) , namely 

(e, V\ M g , ai ,a 2 ,...,an(q,Pl,P2, ■ ■ ■ ,Pn)) - Qstf ^12) (<?) I M aim> ...,a n (Pi, P2, ■ ■ ■ , Pn)) , (116) 

but now the soft current J^\{q) * s 110 longer equal to the eikonal current in Eq. fl83|) . 
In fact, since p\ and p 2 can become collinear, the internal partonic line with momentum 
p l + p 2 in M ai ,a 2 ,...,a n (Pi,P2, ■ ■ -,Pn) is close to the mass shell ( (p l + p 2 ) 2 = s l2 — > 0). Near 
the mass shell, soft-gluon radiation from this internal line leads to soft singularities and it 
cannot be neglected. 




(c) (d) 

Figure 10: Soft-gluon insertion diagrams for the soft-collinear limit. 



The explicit expression of the gluon current J^il) can ^ e derived by working in a 
physical gauge and using the soft-gluon insertion rules described in Sect. |3TT] (Fig. [TUp. We 
find 



(pi +P2) 2 



i=3 

+ (Ti + T 2 



Pi-q (pi+P2 + q) 2 

2(p?+rf) 



T 1 



Pi 

Pi ■ q 



P2 

P2 ■ q 



(pi +P2 + q) 2 



117) 



We discuss the three contributions on the right-hand side in turn. The first contribution 
comes from the usual eikonal insertions on the external parton lines i = 3,...,n (the 
diagrams (a) in Fig. p!0|). 
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The second contribution comes from the soft-gluon emission from the external partons 
Pi and P2 (diagrams (b) and (c) in Fig. |1(]). The factor in the square bracket is the 
usual contribution from the eikonal vertices and propagators of the lines p\ + q and p 2 + 
q. The factor in front of the square bracket has the following origin. In Eq. ( |116| ) we 



have already factorized the tree amplitude M. ai ,a 2 ,...,a n (pi,P2 } ■ ■ ■ ,Pn), which contains the 
propagator factor l/(pi + P2) 2 ■ In diagrams (b) and (c) of Fig. [10|, this propagator is 
instead absent, and it is replaced by the propagator l/(pi + P2 + q) 2 of the internal line 
with momentum p 1 + p 2 + q. Thus, the rescaling propagator factor {jpi +P2) 2 / {pi +P2 + q) 2 
has to be applied to the contribution to the current. 

The third contribution is the eikonal factor due to the soft emission from the nearly 



on-shell internal line Pi + p 2 + q (diagram (d) in Fig. [10]) 



Note that we have neglected diagrams in which p\ and p 2 are not produced by a single 
line with momentum p\ +p 2 . These diagrams are not collinearly singular (see the discussion 
in Sect. EO) in the physical gauge we are working on. 



Note also that, as the eikonal current in Eq. (^4j), the soft current in Eq. ( |117| ) satisfies 
the property q^J^(q) = Y^i=\Ti- The ensuing current conservation, which follows from 
Eq. (|31~D, guarantees the gauge invariance of the factorization formula (|116|) . 

Using Eq. (|116|) we could now perform the collinear limit of the tree-level matrix element 
-Mai,a 2 ,...,a„(Pi, p 2 , ■ ■ ■ , p n ) on the right-hand side. However, since we are eventually inter- 
ested in the soft-collinear limit of the square of the matrix element M.g t ai,...,an(<liPii • • • ,Pn), 
this procedure is not convenient for two reasons. First, we have to introduce collinear split- 
ting functions for the various colour subamplitudes that contribute to the colour vector 
I •M. ai ,a 2 ,...,a n (Pi,P2, ■ ■ ■ ,Pn))- These splitting functions differ from the Altarelli-Parisi split- 
ting functions of Sect. |2.1| (roughly speaking, the former are the square root of the latter) 
and, although they are well known 0], we shall show that they are not really necessary for 
the final result. Secondly, the colour-charge transformation produced by the soft current 
in Eq. (|116| ) implies a non-trivial relation between the colour-subamplitude decomposition 
of the matrix element on the left-hand side and the corresponding decomposition for the 
matrix element on the right-hand side. This non-trivial relation complicates the colour 
structure and leads to mixed soft-collinear splitting functions |2l[| , whose introduction can 
instead be avoided or, at least, simplified. 



In other words, if we square the right-hand side of Eq. ( |116| ), the soft current J 



'(12) 

produces non-trivial colour correlations of the type T 1 ■ Ti or T 2 ■ Tj (with % = 3, . . . ,n) 
between A4(pi,p 2 , . . . ,p n ) and M j( (pi,p 2 , . . . ,p n ). Thus, we cannot perform the collinear 
limit k± — > by simply using the known 0(ao) results of Sect. |2.1| for the colour-summed 
squared amplitude |.MQoi,p2, • • • ,Pn)\ 2 - 

The whole procedure can be simplified by exploiting the QCD coherence properties of 
soft-gluon emission. We rewrite Eq. ( |1 1 7| ) by splitting the soft current in two terms as 
follows: 

JU") = J>. ^ + (T, + T 2 ) + S J Wa) , (118) 
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where 



x TV f \ (Pi + P 2 " 



(Pi +P2 + g) 2 



l 1 r J-2 K 1 1 + J- 1 



Pi-q P2-q (pi+P2)-q 



;ii9) 



The two terms, 5 J (12) an d the other contribution on the right-hand side of Eq. ( |118| ), are 
separately conserved and, thus, the decomposition in Eq. ( |118[ ) does not spoil the gauge 
invariance. 

Then we note that each of the two terms in Eq. ( |118| ) has the correct scaling behaviour 
of 0(l/q) when q — > 0. Their collinear behaviour is nonetheless quite different. Performing 
the limit k± — > at fixed k\jq in Eq. ( |119[ ), the propagator factor {j>\ +P2) 2 / (pi +P2 + q) 2 



is of 0(1) but the term in the square bracket is of 0(k±/q). Thus, the contribution of 
5.7(12) to the soft current J (12) is suppressed by a relative factor of 0(k±) in the collinear 
region, and it can be neglected in the soft-collinear limit. 



We conclude that in the factorization formula we can consistently use the following 
approximation for the soft current in Eq. ( |117f ): 

J fa (g) * E Ti ^- + T (12) V \ -J Pa , (120) 
[ ' 7^3 Pi-V (Pi+P2)-q 

where Th 2 ) = Ti + T 2 - The subdominant effect of S J V12) is due to the cancellation between 
the different contributions in the square bracket on the right-hand side of Eq. (|119|) . The 
cancellation is a typical consequence of colour coherence. When the parton momenta pi 
and p 2 become collinear, they radiate soft gluon in a coherent way, i.e. as a single parton 
with momentum p\ +p 2 and colour charge T( 12 ) = T\ + T 2 (see the last term in Eq. ( |120[) ). 



The expression in Eq. ( p. 20] ) is certainly simpler than that in Eq. ( |117|) . More impor- 
tantly, it depends on the colour charge T(i2) rather than separately on the colour charges 
Ti and T 2 . This implies that, when we square the amplitude in Eq. (|116| ), the partons p\ 
and P2 are no longer colour-correlated, and the collinear limit k± — > can be performed by 
using the collinear factorization formula (0). We obtain the final soft-collinear factorization 
formula: 

2 

\Mg, au a 2 ,...,a n {q,Pl,P2, ■ ■ ■ ,Pn)\ 2 ^ (47T/i 2e a S ) 2 

$12 

■ {M a ,...,an(p, ■ ■ ■ ,Pn)\ Paxaz [ J \l2)tx (?) J (12) (?)] I M a ,...,a n (P> ■ ■ ■ , Pn) ) , (121) 

where the matrix elements on the right-hand side are obtained by removing the soft gluon 
q and by replacing the partons a\ and a 2 by the single parton a that leads to the collinear 
splitting process a — > ai + a2- Since these matrix elements are vectors in the colour+helicity 
space, both spin and colour correlations are present in Eq. ( |121| ). 

The spin correlations are exactly the same as in Eq. ([?]). The spin indices s, s' of the par- 
ent parton a are correlated by the Altarelli-Parisi splitting functions P ai a 2 = Paia 2 ( z > ^-J-! e ) 
in Eqs. 



The colour correlations affect all the partons and are analogous to those in Eq. 
They are produced by the square of the soft current: 
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(12) ((Pi+P2)- q y 



~ T i -T j S ij {q) + 2Y d T i -T {l2) S i{12) {q) . (123) 

i,j=3 i=3 

In Eq. ( |123| ) we have neglected the last term on the right-hand side of Eq. ( |122| ) , because it 
is not collinearly singular. We have also introduced the eikonal functions S,ij{q) of Eq. ( [38|) 
and the analogous eikonal function Sin 2 )(q), 

2(fiii 



«5i(12)(?) 



(124) 



■Sjgl-Siq + S2g) 

From Eqs. ( |121| ) and ( |123| ) we can see that the soft-collinear limit at 0(a\) is simply and 
fully described in terms of the same factors, namely, soft eikonal functions and Altarelli- 
Parisi splitting functions, which control the soft and collinear limits at O(as), respectively. 



The soft-collinear limit at 0(a§) was first studied by Campbell and Glover |^T]|. They 
neglected spin correlations and considered the singular behaviour of the colour subampli- 
tudes. This behaviour, which was extracted by directly performing the singular limit of 
known squared matrix elements, was given in terms of two different factors. The first factor 
(see Sect. 4.4 in Ref. |2l|) refers to subamplitudes in which the collinear partons are not 
colour-connected and it corresponds exactly to the 5^ (g) -term in Eq. (|123|) . The second fac- 
tor regards the subamplitudes in which the collinear partons p\ and p 2 are colour-connected. 



This factor is given in Sect. 5.2 of Ref. pi] and it can be written as 



S; 



2(sii + s i2 ) 



i;ql2 



>iq 



Sin 



Z + 



Si q + ZS\ 2 
S\2q 



2(sa + s i2 ) 



"iq 

4(fq + s i2 ) 

Siq(Slq + S 2q ) 



$lq + S 2g 



1 + 



jU2_ 

S l2q, 



Z 
i s lq 



$lq + S 2q 



(125) 
(126) 
(127) 



Since z is the longitudinal momentum fraction carried by pi in the collinear region, the term 
in the square bracket of Eq. ( |126j ) vanishes in the collinear limit and Eq. ( |127| ) follows. This 
simplification, which is due to colour coherence, was not performed in Ref. [^TJ. Taking 
it into account, we have S i]q i 2 ~ 2iSj( 12 )(g), which, when inserted in Eq. ( |123| ), shows the 
equivalence of our results with those of Ref. |2l 



Our derivation of the soft-collinear factorization formula ( |121| ) can straightforwardly be 
extended to higher orders. We can consider the limit where a single gluon with momentum 
q becomes soft and, at the same time, m partons, say pi, ■ ■ ■ ,p m , become simultaneously 



collinear (see Sect. |2.3|) . In this limit the factorization formula is 



\Mg, ai , 



Mr, 



, an (q,Pi, 

■ ,a n {xp,.. 



,Pr, 



,Pn) 



m— 1 



t 

(l...m)ju 



(l) J (l...m) 



Si.. 



M a ,...,a n (xp,..., Pn )) ,(128) 
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where P ai ...a m = Pat.. am 1S the spin-dependent splitting function in Eq. (p5|) , and the soft 
current J n... m )(q) is: 



J?i... m) (<i) - E r< — + r ( i-m ^ 1 + '-- p " , (129) 



i=m-\ 



with T(i... m ) = Ti + . . . + T m . Squaring the soft current as in Eqs. ( |122j ) and ( |123| ), we 
obtain 

n n 

j\i... m ),(Q) J (l...rn)(Q)- E T i -T j S ij {q) + 2 E P t -T {1 ... m) S l{1 ... m) (q) , (130) 

i,j=m+l j=m+l 

where <S i( i... m )(g) = 2(s a + . . . + s <m )/[s i9 (si, + . . . + s mq )). 

The proof of these results is very simple. The soft current J(i... m )(q) is derived by 
using the soft-gluon insertion rules as in Eq. ( |117| ). Then, the coherence argument used in 
Eqs. (|118|) and ( |119| ) can iteratively be applied to any vertex in the m-parton dispersive 
amplitude V ai ... am of Eq. p7|). This leads to the expression in Eq. ( pL29| ) . 



3.5 Multiple soft and soft— collinear limits 



In Sects. (|3.1|) and ( |3.3j ) we have discussed in detail single and double soft-gluon emis- 
sion. The generalization to multiple soft-gluon radiation is straightforward. If we consider 



the matrix element M., 



,9,0-1, 



, . . . ,p n ) when the k gluons with momenta 



qi,- 

Eq. 



qt become soft simultaneously, we can still write a factorization formula similar to 
by performing the simple replacement 



2 2e 

9s^ 



ja\. 



Qk) 



(131) 



where ax, . . . , and \i 



. , f/,k denote the colour and Lorentz indices of the soft gluons. 
As the two-gluon current ^^(91,52) i n Eq. ( |101| ), the multigluon current on the right- 



hand side of Eq. (|131|) is obtained by working in a physical gauge and using the soft-gluon 
insertion rules described in Sect. |3~T| . Of course, the explicit expression of J(qx, ■ ■ ■ ,qk) 
turns out to be quite involved, because all the possible interactions between the soft gluons 
have to be included without using any soft approximation. 

It is also clear that the soft-gluon insertion rules can be used to derive a factorization 
formula analogous to Eq. ( |116|) for the multiple soft-collinear limit in which k gluons are 
soft and m partons become collinear simultaneously. More importantly, it is worth while 
pointing out that the coherence argument leading to Eqs. ( |121| ) and ( |128| ) still applies. 
Thus, the factorization formula can be written as 



\M 9 ,.... 



g,a\,...,a m , 



,a n (qU---,qk,PU---,Pm,---,Pn)\ 2 - (~ ^TT ^ a S ) ' 



Sl...m 



m—1 



(132) 



• (M (xp, ..,p n ) I P 0,1. ..dm (l...m) a,...,a n 

(xp, ..,p n )) 



Equation (|132|) does not involve any additional factor with respect to those that are nec- 
essary to deal with the multiple collinear and multiple soft limits separately. The spin- 
dependent splitting function P ai ...a m = Pat. .a * s exactly the same as that in Eq. 
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The current Jn m )(qi, ■ ■ ■ , qk) is completely analogous to the soft current J(qi, . . . , qk) in 
Eq. ( |131| ). As the latter, the former is constructed by inserting the soft gluons only on the 
(1 +n — m) external parton lines with momenta p 1 + . . . + p m , p m +i, ■ ■ ■ , p n , an d each inser- 
tion on the collinear parton p\ + . . . + p m is taken into account by the simple eikonal factor 
(Ti + . . . + T m )( Pl + . . .+PmY/(pi + ■ ■ ■+p m ) ■ q, despite the fact that ( Pl + . . . +p m f ^ 
(see e.g. Eqs. ( |129| )). 

In the most general case, the infrared singularities of the tree-level amplitudes are 
produced by the multiple collinear decay of hard partons and by the associated radiation 
of soft gluons and gg-pairs. The corresponding factorization formula can be constructed in 
a straightforward manner by using the rules derived and illustrated throughout the paper. 
Using a shorthand symbolic notation, we have (Fig. [11]) 

\M\ 2 ~ (M hard \ (liP^I S\M hard ) . (133) 

Here Aihard denotes the factorized amplitude that depends only on the momenta of the 
hard partons. The factor flj Pi — IL P% 1 1 is the product of the spin-dependent splitting 
functions for the collinear decay of i = 1, . . . , I hard partons. The factor S is a colour 
matrix that takes into account the radiation of soft partons. It has to be computed exactly 
at the tree level, but its external gluon lines are coupled to the hard partons by using the 
eikonal approximation as in the calculation of the current J(i... m )(?i, • • • , qk) in Eq. (|132p . 
Note that spin and colour correlations are factorized independently. This decoupling follows 
from colour coherence. 




Figure 11: General structure of infrared factorization at any perturbative order. 
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4 Summary 



In this paper we have studied the infrared structure of tree-level QCD amplitudes in all 
the possible soft and collinear limits. 

We have first considered the collinear behaviour. We have shown that, in the limit in 
which m partons become parallel simultaneously, the singularities are given by the universal 



factorization formula fl25|) and are controlled by process-independent splitting functions 
that generalize the customary Altarelli-Parisi splitting functions. These splitting functions 
fully take into account the azimuthal correlations produced in the collinear decay. We have 
presented a recipe to compute the splitting functions at any perturbative order and we 
have performed their explicit calculation at 0(a§). 

Then we have studied the soft behaviour and shown how to construct soft factorization 
formulae at any order in ag- We have considered the limit in which a qq pair becomes 
soft and we have computed the corresponding singularity at 0(a\) in terms of a simple 
universal insertion factor. We have then recalled the known results about the limit in which 
two gluons become soft. This limit is controlled by an 0(a§) soft current that is tensor in 
colour space and generalizes the eikonal current at 0(a>s)- We have obtained a compact 
expression for the square of the two-gluon current that, in particular, shows the absence of 
colour correlations in the case of four- and five-parton amplitudes. 

Finally, we have studied the mixed soft-collinear limit and pointed out that its descrip- 
tion does not require the introduction of new infrared factors. Exploiting the coherence 
property of soft gluon radiation, we have been able to show that the singularities are given 
by a factorization formula written only in terms of the soft currents and of the splitting 
functions that control the soft and collinear limits, respectively. 

These results are one of the necessary ingredients to extend QCD predictions at higher 
perturbative orders. In particular, our calculation of the 0(a\) singular factors is relevant 
to setting up general methods to compute QCD jet cross sections at NNLO. 
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Appendix: Soft limits of four- and five-parton ampli- 
tudes 



In general, soft factorization formulae involve colour correlations. As shown in Eqs. 
( |112|) and ( |123|) , at 0(a§) the correlations are completely given in terms of products of 



colour-charge factors Tj • Tj. 

In this appendix we collect the factorization formulae for the 0(a\) soft (and soft- 
collinear) limits of the square of the matrix elements with four and five partons plus an 
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arbitrary number of colourless particles. In these particular cases, using colour conservation, 
it is possible (see e.g. the Appendix A in Ref. 0] ) to express the products Tj • Tj in terms 
of the Casimir invariants Cj (Cj = Cp if i = q, q and Cj = Ca if i = g) of the hard partons. 
Thus the colour algebra completely factorizes and colour correlations cancel. 

Note that two of the hard partons in the four- and five-parton amplitudes necessarily 
form a particle-antiparticle pair a, a. This further simplifies the combinations of Casimir 
invariants that appear in the factorization formulae. 

Emission of a soft gg-pair 

We consider the four-parton amplitude M. q> q, a ,a{qi-i Q2,Pi,P2) in the limit in which 



<Zi><72 — ► 0. From Eq. (|95Q , we get 



\M q ,g, ara (qi,q 2 ,Pi,P2)r - {^fi 2e a s YT R C a Xu+X 22 -2X 12 \M a - a { P i,p 2 )\ 2 • (A.i; 



In the case of five partons we get 



\M q ^a,a 3 (?1 , 92, Pi , Pi, Ps) 



Ca {111 + ^22 — 2Xi2) 



(A.2) 



+ Ca 3 ( X33 + X12 — X13 — X 



23; 



\M {Pl,P2,Ps)\ ■ 



4 C a S 12 ( qi ) S 12 (q 2 ) + C A (2 «Si 2 - <Sn - S 22 

(A.3) 



The soft function Xy = %ij(qi, q 2 ) is given in Eq. (j96j). 

Emission of two soft gluons 

We consider the amplitude M. g>g>aj a (qi, q 2 ,Pi,P 2 ) in the limit in which the two gluons 
become soft. Using Eq. ( |112j ), we get 

\Mg ig ^ a (qi,q 2 ,Pi,p 2 )\ 2 ^ (47r/i 2e a s ) 2 C a 
In the case of five partons we get 

\M g ^ a> a,a 3 (<7l,<?2,Pl,P2,P3)| 2 - (47T/i 2e a S ) 2 \M afi ,a 3 (Pi , P2 , P3) 1 2 

(2C a - C a3 ) S 12 ( qi ) + C 03 (5 13 (gi) + <S 23 (gi)) 

(2C a - C a3 ) S 12 (q 2 ) + C a , («S 13 (g 2 ) + S 23 (q 2 )) (A.4) 



Ca (2<5>12 — Six — S 22 ) + C a3 (<Sl3 + 1S23 — S33 ~ S 



'12/ 



The soft functions Sij(q) and <Sy = <Sy(gi, g 2 ) are given in Eqs. flSS|) and ( |109| ), respectively. 
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Soft— collinear limit 

We consider the amplitude M.g^ ai ,a 2 ,a 3 (q,Pi,P2,P3) i n the limit in which q — > and 
s 12 — > 0. Using Eq. (|121|) , we get 

\M g , aua2ta3 (q, Pl , P2 ,p 3 )\ 2 ~ — (4vr/i 2e « s ) 2 C a3 S 3(12) (q) P*f a2 T a s a s ^xp,p 3 ) . (A.5) 

•512 

In the case of five partons we get 

\M g , aua2 , a3 , a4 (q 1 ,p 1 ,p2,P3,PA)\ 2 - — (4vr/i 2e a s ) 2 T a s a s ' {xp,p 3 ,p 4 ) ^ (A.6) 
• (C as + C a4 - C a ) S 34 (q) + (C a + C a3 - C ai ) S 3(12) (q) + (C a + C a4 - C a3 ) S 4(12) (q) . 

Here a denotes the parton that decays collinearly, a — > 0102, T£ s (xp,...) is the spin- 
polarization tensor in Eq. (0) and i^fao ^ s ^ ne spin-dependent splitting function in Eq. (0). 
The soft functions Sij(q) and «Si(i2) (?) are given in Eqs. ( |8^ ) and ( |124| ), respectively. 
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